Sideward flow in Au + Au collisions between 2A and 8^4 GeV 
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Abstract 

Using the large acceptance Time Projection Chamber of experiment E895 at 
Brookhaven, measurements of collective sideward flow in Au + Au collisions 
at beam energies of 2, 4, 6 and 8^4 GeV are presented in the form of in- 
plane transverse momentum (p x ) and the first Fourier coefficient of azimuthal 
anisotropy v\. These measurements indicate a smooth variation of sideward 
flow as a function of beam energy. The data are compared with four nu- 
clear transport models which have an orientation towards this energy range. 
All four exhibit some qualitative trends similar to those found in the data, 
although none shows a consistent pattern of agreement within experimental 
uncertainties. 

PACS numbers: 25.75.Ld 
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Sideward flow was the first type of collective motion to be identified among fragments 
from relativistic nuclear collisions Q. It consists of a preferential emission in the plane 
defined by the incident nuclei (the reaction plane); at relativistic energies, nucleon emission 
towards the projectile side is favored forward of the center of mass rapidity, while the target 
side is favored at backward rapidities. This behavior is normally attributed to a release of 
compressional energy, and thus is sensitive to the integrated effect of the nuclear pressure 
generated in the collision. Models indicate that sideward flow is established during the early, 
high density stage of the heavy ion collision, and that it is minimally distorted during the 
subsequent evolution. 

A Quark Gluon Plasma (QGP) might be formed in heavy ion collisions at sufficiently 
high energies, in contrast to the purely hadronic matter that exists throughout the collision 
process at lower energies. Near the transition between these two regimes, it is argued 
that the increased entropy density leads to a "softest point" in the nuclear equation of 
state (EOS), which could generate a minimum in the pressure-driven sideward flow at the 
relevant beam energy 0-fH- Earlier, it had been suggested that sideward flow at CERN 
energy would have a magnitude that depends on whether or not a plasma is produced || . 
Initial calculations incorporating a softest point in the framework of one-fluid relativistic 
hydrodynamics featured a prominent minimum in the sideward flow for near-central Au + 
Au collisions around 5 A GeV ||. However, the softening effect is reduced when allowance 
is made for the finite size of the hydrodynamic system |J, and it occurs at higher beam 
energies in a three- fluid model |4]]. Most recently, it has been argued that the shape of the 
rapidity dependence of sideward flow may be a QGP signature 0. 

In the past, hydrodynamic calculations have frequently been the first to predict new 
collective phenomena, whereas microscopic transport models have typically reproduced flow 
measurements subsequently with better agreement. Relativistic transport codes yield mea- 
surably different flow amplitude near the "softest point" beam energy, depending on whether 
or not a phase transition is simulated [BUD!. Furthermore, the smooth excitation function 



for elliptic flow observed in E895 |L0| has been interpreted as a possible phase transition 
signature because a transport model comparison is consistent with a progressive softening 
of the EOS with increasing beam energy ||. However, a transition from hadronic to string 
degrees of freedom has since been put forward as an alternative interpretation for such an 



EOS softening |TT|. There are several further reasons why the 2-8 A GeV energy range is 
especially interesting: it is largely unexplored; many inelastic NN channels open up within 
this relatively narrow range of beam energies; and models suggest that the highest baryon 
density is reached in this region [p~2]j . 

We report proton sideward flow measurements for Au + Au collisions at kinetic energies 
of 1.85, 3.9, 5.9, and 7.9A GeV in experiment E895 |13[ at Brookhaven's AGS. The data pre- 
sented come from the main E895 subsystem — the EOS Time Projection Chamber (TPC). 
E895 allows a seamless extension to higher beam energies of the detailed flow excitation 
functions already measured [|14| , |T5"| using the same TPC at the Bevalac. The TPC offers 
good acceptance for charged particles over a substantial fraction of Atc solid angle, as well 
as particle identification via energy loss measurement, although ambiguities in separating 
protons from both pions and deuterons increase with beam energy. The full event recon- 
struction capability of E895 allows determination of flow even when the correlations used to 
establish the reaction plane orientation are relatively weak. 
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The estimated reaction plane azimuth $ for an event is based on the orientation of 
Q = J2j w pf/pfy where j runs over all baryonic fragments in the event, is momentum 
in the plane perpendicular to the projectile direction, and we use the weighting factor w = 
yj/max(|y'-|, 0.8), where = ?/j ab /y mid — 1 is the center-of-mass frame rapidity for fragment 
j, and y mid is half the rapidity gap between target and projectile JjiJ. Thus, y' denotes 
normalized rapidity such that the target and projectile are always at y' — — 1 and +1, 
respectively. Detector asymmetries and inefficiencies result in anisotropies in the determined 
reaction plane azimuth $. Uniform distributions in $ are recovered after each track is 
assigned a weight according to its y and magnitude and direction. 

The centrality of collisions is characterized in terms of charged ejectile multiplicity M as 
a fraction of M max , the value near the upper limit of the M spectrum where the height of the 
distribution has fallen to half its plateau value |T7]]. The events used in our analysis come 
from the region where sideward flow is at or near its maximum — multiplicities between 0.5 
and 0.75 times M max . Models indicate that this region corresponds to impact parameters 
distributed mostly between 5 and 7 fm. After centrality selection, the analyzed samples 
contain 12, 24, 9.3 and 7.2 thousand events at 2, 4, 6 and 8A GeV, respectively. 

The mean proton transverse momentum projected onto the reaction plane, (p x ), is pre- 
sented as a function of rapidity in Fig. [I]. Using the prescription described in Ref. (JTS 



these and all subsequent flow signals are corrected for finite resolution in determining the 
reaction plane. Dispersion correction factors are 0.89, 0.79, 0.62 and 0.43 at 2, 4, 6 and 
8A GeV, respectively. In E895, there are known distortions below p± ~ 0.3 GeV/c caused 
by track reconstruction inefficiencies and by the breakdown of proton - 7r + separation at 
some rapidities. At Bevalac beam energies, where these distortions are not a factor, we 
have studied the contours on scatter plots of the components p x and p y within various 
rapidity gates, after events are rotated so that estimated reaction planes are aligned with 
the x axis. It is observed that the contours are concentric, and (p x ) is constant in any slice 
of p y , with deviations of less than 5%. Therefore, full acceptance in the p± plane is not re- 
quired to extract the sideward flow, provided an appropriate p y cut is applied to remove the 
problematic region. For p y > 0.3 GeV/c, the expected flat behavior in (p x (p y )) is observed 
at E895 beam energies. We assign these plateau values to (p x ), with the assumption that 
the behavior described above is a general property of sideward flow which does not change 
between Bevalac and AGS energies. Moreover, GEANT-based simulations of the detector 
response without this assumption yield corrected {p x ) results that are consistent within the 
reported uncertainties. 

Shapes of (p x {y)) are normally close to linear over an interval centered on midrapidity, 
and a function Fy' + Cy' 3 typically fits the {p x (y')) distribution over the y' region dominated 
by participant fragments. It has become common to use the fitted linear coefficient F (or 
F y = F/y mid , the corresponding slope for unnormalized rapidity) to characterize the overall 
strength of the sideward flow. We average the fitted coefficient F with and without imposing 
C = 0, and the difference generally dominates the systematic uncertainty in the slope (which 
is large compared with the statistical error). Figure ^] presents both F and F y as functions 
of beam energy, along with the same quantities for comparable centrality, as measured in 
the same detector at lower energy fC5" |, and in E877 (19| at maximum AGS energy. Nucleon 
rapidity spectra dN/dy' at fixed centrality have close to the same form at different beam 
energies, and this observation motivated the use of normalized rapidity y' in the definition 
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of F. Both flow observables decrease steadily with increasing beam energy over the E895 
range. The significant decrease in F y at E895 energies (in contrast to a flat or slightly 
decreasing trend in F y (E] 3eam ) at Bevalac energies fll5|), is interpreted in the hydrodynamic 
picture as an increased deviation from ideal fluid behavior (constant F y (E heam ) pO|Jl|), most 
plausibly viscosity arising from the increasingly copious particle production in this region 
f20fl . Phenomenologically, the observed trends in F and F y between 0.2 and 10v4 GeV can 
be related to the steady decrease in the azimuthal anisotropy of the proton distribution v± 
(defined below) over this range, in combination with the variation of (p±), which increases 
steeply at Bevalac/SIS energies but appears to approach saturation at higher beam energies 

E- 

Within a rapidity window, flow causes anisotropic distributions of track azimuths <fi 
relative to the reaction plane. These anisotropies generally can be well described by the 
truncated Fourier expansion 

dN/d(j) « v [1 + 2vi cos <\> + 2v 2 cos 20] . (1) 

The first Fourier coefficient, v\, reflects the azimuthal angular part of the sideward flow 
correlation, and is related to the (p x ) sideward flow observable according to 

(p x ) = ^ J v 1 (p ± ) P ^dp\ (2) 

The V2 coefficient in Eq. (|Ij) represents elliptic flow |I] , already reported for the present E895 



data set |LtJ . Figure y presents measured V\ coefficients for protons as a function of rapidity, 
at the four E895 beam energies. p± gates, as labeled, were applied when generating the 
vx(y) spectra reported in Fig. [| (but p± gates were not used in any of the (p x ) analyses). 

The data in Figs. [I] through |3| do not show evidence of a dip in the flow excitation function 
H|6|||[|], and previous measurements at both higher and lower beam energy are consistent 
with a smooth extrapolation of the E895 data. The (p x {y)) slopes show no evidence of 
a decrease near midrapidity, while a small flattening effect [f2~l~1 ,|7f appears in Vi(y) at the 
highest energies. The decrease in V\{y) slope at midrapidity is very prominent at CERN 
energy 



22 



Nuclear transport models describe the nuclear collision in part as successive point-like 
nucleon-nucleon interactions. The NN cross sections are mostly taken from experiment and 
include inelastic processes, e.g., production of resonances, pions, etc. These models have 
been successful in reproducing a large fraction of the published flow measurements to within 
a few tens of percent or better. The representative transport models now available include 
RQMD J23[ (Relativistic Quantum Molecular Dynamics), UrQMD H (Ultrarelativistic 



Quantum Molecular Dynamics), ART || (A Relativistic Transport model) and BEM || 
(Boltzmann-Equation Model). In one operating mode of RQMD ("cascade"), the hadrons 
and resonances propagate freely between binary collisions, and the equilibrium pressure is 
close to that of an ideal gas. RQMD also contains an option (the "mean-field" mode) which 
allows additional pressure to be generated in the high density stage. UrQMD provides only 
a cascade mode at these energies, but features a completely independent implementation 
from RQMD. BEM is based on the relativistic Landau theory of quasiparticles. In addition 
to a cascade mode, two types of momentum-dependent EOS can be selected in BEM: a 
scalar potential (soft EOS, K = 210 MeV) and a vector potential (hard EOS, K = 380 
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MeV). In the ART model, there is likewise provision for a soft and stiff EOS, but using 
different phenomenological prescriptions from those in BEM. An early version of the BEM 
code was used in refs [0,0 ; the version reported here yields less sideward flow than before 
25fl , but almost the same elliptic flow [10 125|. All comparisons presented here are subject 



to a systematic uncertainty arising from the fact that all of the model calculations neglect 
formation of composite nuclear fragments, while the flow measurements are for free protons 
only. 

Although a relatively complex pattern of disagreement is observed between data and 
all four models, several general conclusions are suggested by the comparisons. The cascade 
modes of BEM, RQMD and UrQMD all exhibit less proton flow F than observed. The 
significant differences among the codes in cascade mode (up to a factor of 2 or more) indicates 
that the binary scattering part of transport simulations, usually considered to be better 
understood than the "long-range" part, remains a source of significant uncertainty in model 
calculations at these energies. The soft EOS in ART and BEM, and RQMD's mean field, 
all come close to reproducing the F measurements. The relative variation of v\ for data and 



all models as a function of increasing beam energy is suggestive of a softening trend pO 
However, there is a marked tendency for the v i data to favor substantially softer equations 
of state than F measurements. The coefficient v i reflects only the azimuthal angular part 
of the sideward flow correlation, while (p x ) and F also include the effect of p± magnitude 
correlations; thus, the tendency noted above indicates that the models consistently have too 
small a p± magnitude flow correlation relative to their azimuthal angle flow correlation. It is 
also noteworthy that no mode of any of these models is close to simultaneously reproducing 
the E895 elliptic flow [PH , P5|| and our data for F and v±. 

In summary, we report measurements of sideward flow in collisions of heavy nuclei in the 
previously unexplored region between maximum Bevalac/SIS energy and maximum AGS 
energy. Sideward flow decreases smoothly over the 2A to 8A GeV range, and extrapolations 
are consistent with existing measurements at both lower and higher beam energies. A 
new trend of decreasing flow sets-in near the low end of the studied beam energy range 
- the F y excitation function changes from flat or slowly decreasing to a steeper rate of 
decrease, while the slope of the F excitation function changes sign. This change roughly 
coincides with the onset of copious particle production, and is reproduced qualitatively 
by some transport calculations. As at other energies [|IJ, transport simulations in cascade 
mode consistently yield less flow than is observed. The model flow calculations for E895 
energies are subject to significant systematic uncertainties, as inferred from variations among 
different models and from the overall extent and pattern of agreement with experiment. 
At present, these uncertainties appear to be larger than the magnitude of flow signatures 
of physics importance, and so strongly motivate further transport model investigation and 
development focused specifically on this unique energy domain where many new inelastic NN 
channels open up, and meson yields increase steeply. Overall, our flow excitation function 
measurements offer important constraints on possible conclusions regarding a QGP phase 
transition at or above E895 energies. 

We thank the authors of refs. [3,|9 p3l , [2"4f for providing model codes and/or calculations. 
We also thank P. Danielewicz, A. Poskanzer and S. Voloshin for valuable input and discus- 
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FIG. 1. Average proton function of normalized rapidity, y' (the target and projectile 

are always at y 1 = — 1 and +1, respectively). The closed symbols are direct measurements and the 
open symbols are generated by reflection about midrapidity. Histograms are RQMD calculations 
in cascade mode (dashed line) and mean field mode (solid line). 
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FIG. 2. Proton flow magnitude as a function of beam energy; the lower right panel shows 
the measured F y , while the other three panels show identical measurements of the parameter 
F, with different transport model calculations superimposed. The error bars include systematic 
uncertainties. 
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FIG. 3. Fourier coefficients v\ as a function of normalized rapidity, y'. The closed symbols 
are direct measurements and the open symbols are generated by reflection about midrapidity. 
The labeled transverse momentum gates are in units of GeV/c. Histograms are transport model 
calculations, as labeled at the top of each column, and each column of four panels shows identical 
measurements. 
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